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Large-scale patterns, which are well-known from the spiral defect chaos regime of ther¬ 
mal convection at Rayleigh numbers Ra < 10^, continue to exist in three-dimensional 
numerical simulations of turbulent Rayleigh-Benard convection in extended cylindrical 
cells with an aspect ratio F = 50 and Ra > 10^. They are uncovered when the turbulent 
fields are averaged in time and turbulent fluctuations are thus removed. We apply the 
Boussinesq closure to estimate turbulent viscosities and diffusivities, respectively. The 
resulting turbulent Rayleigh number Ra*, that describes the convection of the mean 
patterns, is indeed in the spiral defect chaos range. The turbulent Prandtl numbers are 
smaller than one with 0.2 < Pr* < 0.4 for Prandtl numbers 0.7 < Pr < 10. Finally, we 
demonstrate that these mean flow patterns are robust to an additional finite-amplitude 
side wall-forcing when the level of turbulent fluctuations in the flow is sufficiently high. 


1. Introduction 

The formation of regular patterns close to the onset of a hydrodynamic instability 
in spatially extended flows is well documented for generic cases. The most prominent 


examples are convection rolls in Rayleigh-Benard flow (Busse (1978), Bodenschatz et al. 


(2000)) heated from below and cooled from above, Taylor vortices in Taylor-Couette flow 
(Andereck et al. (1986)) between two rotating concentric cylinders and inclined turbulent 
stripe patterns in plane-shear flows driven by a pressure gradient or a wall movement 


(Barkley & Tuckerman (2005), Duguet & Schlatter (2013)). Specifically linearly unstable 
systems, such as Rayleigh-Benard convection, with a sharp transition threshold to the 
convective flow state allow then for a perturbative expansion about the first unstable 
mode at onset. The expansion leads to an amplitude equation which is simpler than 
the original fluid equations and describes the formation of simple patterns as a function 


of the system parameters (Cross & Hohenberg (1993), Hoyle (2006)). The derivation of 


nonlinear phase diffusion equations allows to model increasingly complex patterns, such 


as spirals or defects (^ 

doyle 

(2006 

)), which have also been detected in experiments, e.g., by 

Croquette 

(1989 

) anc 

Morris et al. ( 

1991). Defects are imperfections in the patterns such 

as dislocations ( 

Cross & Greenside 

(2009)). In large-aspect ratio cells, this results in a 


state of slowly evolving spirals and defects which is known as spiral defect chaos (SDC). In 
such SDC regimes, all symmetries of the governing equations have been spontaneously 


broken (Busse (2003)). For example, the azimuthal symmetry of the roll patterns is 
broken in an extended cylindrical cell. When the temperature difference across the fluid 
layer is further increased, the fluid motion crosses over from the weakly nonlinear to the 
turbulence regime. 

The dimensionless Rayleigh number Ra = gaATH^ /(vk) describes the thermal driv¬ 
ing. It contains the acceleration due to gravity, g, the thermal expansion coefficient at 
constant pressure, a, the outer sustained temperature difference across the layer, AT, 
and the height of the layer (or convection cell), H. The kinematic viscosity of the fluid v 
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and the thermal diffusivity k of the temperature field form the second important dimen¬ 
sionless parameter which relates dissipation in the working fluid to thermal diffusion, the 
Prandtl number Pr = vj k. Turbulence is characterized by an irregular, stochastic and 
three-dimensional fluid motion. Does this however imply that the patterns for velocity 
and temperature which are documented in the SDC regime at lower Rayleigh number 
disappear? If not, how can these patterns be extracted? How robust are they with respect 
to variations of the Prandtl number? And finally, how robust are they with respect to an 
additional side wall-forcing which is added to the momentum equation? These are the 
questions which we want to address in the present work. 

Our investigation is based on three-dimensional direct numerical simulations (DNS) in 
very large aspect ratio cells which are comparable to laboratory experiments in pattern 


formation (Bodenschatz et al. (20001). We conducted a series of DNS of the Boussinesq 
equations in cylindrical cells with an aspect ratio of P = D/H = 50 with D being the 
diameter of the cell. Our investigation extends previous numerical studies of turbulent 


Rayleigh-Benard convection by Hartlep et al. (2005), Hardenberg et al. (20081 and Bailon- 


Cuba et al. (2010) to very large aspect ratios P (see also Chilla & Schumacher (2012)). We 


show that patterns very similar to SDC continue to exist into the soft turbulence regime, 
up to Rayleigh numbers Ra = 500 000 which were accessible here. The patterns are 
covered by an increasing amplitude of velocity fluctuations and become visible only after a 
time averaging over a sequence of flow snapshots. We then determine turbulent viscosities, 
and diffusivities, /c*, for the time-averaged roll patterns and define turbulent Rayleigh 
and Prandtl numbers for the mean flow patterns. That means we replace 

Ra —)■ i?a* = - -<C Ra and Pr —)■ Pr* = — ~ 0(1) . (1-1) 

These turbulent Rayleigh numbers are significantly smaller than the original ones and 
fall consistently back into a range that corresponds with the original spiral defect chaos 
regime. The turbulent Prandtl numbers decrease as well and remain smaller than one. 
The value of Pr* depends weakly on Pr. It increases with increasing Pr. 


2. Numerical model 

We perform DNS of the three-dimensional Boussinesq equations which are given by 


du^ 

dt 


dT 

l)t 


du, _ 1 dp d^Ui (rr ^ f 

Jb J rU 

^=0 

dx, ’ 

dT 


d'^T 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


where p(xk,t) is the pressure, Ui(xk,t) the velocity field, po the constant mass density, 
and T{xk, t) the temperature field. The temperature Tq is a reference temperature. Sum¬ 
mation is applied over index j in Eqns. (2.1) - (2.3), where i,j,k = 1,2,3. The last 


term on the right hand side of Eq. (2.1), which is denoted by fi{xk,t), stands for an 


additional volume forcing. It is applied close to the side walls of the convection cell and 
designed such that it enforces the azimuthal symmetry in the vicinity of the side walls. 
This additional forcing is applied for two runs only which are discussed in section 3.3. 

The velocity field has a no-slip boundary condition on all walls. The temperature 
boundary condition is isothermal at the top and bottom plates and adiabatic at the 
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Run 

Ncf, X Nr X Nz 

Ra 

r 

Pr 

Nu 

Re 

Nbl 

Rat 

Pr* 

1 

601 X 401 X 97 

5 000 

50 

0.7 

1.84 

14 

30 

— 

— 

2 

601 X 401 X 97 

500000 

50 

0.7 

7.72 

204 

10 

4500 

0.21 

2a 

1201 X 601 X 141 

500000 

50 

0.7 

7.25 

209 

16 

4800 

0.21 

3 

601 X 401 X 97 

500000 

50 

3 

8.38 

62 

10 

7300 

0.38 

4 

601 X 401 X 97 

5 000 

50 

10 

2.00 

1 

29 

- 

- 

5 

601 X 401 X 97 

500000 

50 

10 

8.71 

22 

9 

39000 

0.40 


Table 1. Parameters of the different simulations. The column Nbl displays the number of grid 
points inside the thermal boundary layer. The last two columns display the turbulent Rayleigh 
and Prandtl numbers, respectively, analyzed for selected runs. Run 2a is conducted at a higher 
grid resolution as run 2. 


side wall. The problem is formulated in cylindrical coordinates (r, (p, z) and solved by a 
second-order finite difference scheme (Verzicco & Camussi ( 2003| )). In Table we list 
the simulation parameters and grid resolutions for each DNS run. The azimuthal spacing 
(A^) is uniform, the radial (A^) and axial (A^) grid sizes are nonuniform. The radial 
mesh gets finer towards the side wall. It is obtained by a geometric scaling relation 
which clusters the grid points less than Tchebychev collocation points when the side 
wall is approached, r —> D/2. This also means that azimuthal grid spacing rA^ grows 
for growing r. The grid resolution in the cylindrical cell with very large aspect ratio 
puts a challenge to the simulations, in particular in terms of the azimuthal resolution. 
Following Grotzbach (19831, we tested our DNS grid by calculating the global maximum 
of the geometric mean, A = max{^/rA//j/A/PA/). The criterion states that A < tttjk, 
whe re vk is the Kolmogoro v dissipation length (see also Emran fc Schm^cher (20081 


and 


Shishkina et al. 


(2010)). The ratio A/px is 2.7 for run 2 in Table and 1.7 for 
comparison run 2a at a higher resolution. Turbulent heat and momentum transfer are 
compared in Tab. 

All runs start with the diffusive equilibrium state which is perturbed randomly. Length 
scales are normalized in units of H, velocities in units of the free-fall velocity Uf = 
gaATH, time in units of the free-fall time Tf = H/Uf and temperatures in units 
of AT. The turbulent heat and momentum transport are measured by the Nusselt and 
Reynolds numbers, respectively. They are given by 


Nu=l + \/ RaPr{u^f)v,t , Re = ■ (2.4) 

The notion {pv.t stands for an ensemble average taken as a volume-time average in DNS 
case. From here on, we will omit the tilde for dimensionless quantities. 


3. Results 


3.1. Observations 

Figure shows a sequence of three-dimensional streamline plots viewed from the top 
for Pr = 0.7 and F = 50. Panel (a) and a magnification in (b) are for Ra = 5 000. For 
this Rayleigh number value almost no difference was found between an instantaneous 
snapshot and the time average which is taken over 100 Tf and not shown in the figure. 
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Figure 1 . Streamlines of the velocity field (view from the top) and contours of the temperature 
field in a Rayleigh-Benard convection cell, (a) Instantaneous velocity field pattern at Ra — 5 000 
and (b) magnification, both taken from run 1 in table (c) Corresponding temperature field 
in mid plane. (d,e) Instantaneous streamline plot and its magnification at Ra = 500 000. (f) 
Corresponding temperatnre field in mid plane. (g,h) Streamline plot of the time-averaged velocity 
field at Ra = 500 000 and magnification, (i) Corresponding time-averaged temperature field in 
mid plane. The time average in (g)-(i) is taken over r = 200T/. Panels (d)-(i) are for run 2 
from tableAll data are for Pr = 0.7. The yellow boxes in panels (b) and (h) highlight defects 
in the patterns. 


The corresponding temperature pattern is displayed in panel (c). Panel (d) of the same 
figure and its magnification (e) display an instantaneous streamline plot at Ra = 500 000. 
Both figures reflect the large amplitude of turbulent fluctuations. The fluctuating nature 
of the temperature field is also obvious in Fig. m- The snapshots appear at a first 
glance almost featureless. The bottom panels (g)-(i) show the time averages, which are 
obtained for a duration of 200Tj, and its magnification. The temperature plots (f) and (i) 
recapture patterns which have been discussed in Hartlep et al. ( |2005 ) for similar Rayleigh 
and Prandtl numbers in rectangular slabs with P = 10. This holds particularly in the 
center of the convection cell. The magnified view of panel (g) in (h) confirms the well- 
known result that the mean flow rolls end perpendicular to the side wall which underlines 
that the grid resolution is sufficient. 

The time-averaged plots (g)-(i) recapture now patterns that are similar to SDC, i.e., 
to those which are observed in panels (a)-(c) of the same figure for the Rayleigh number 
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Figure 2. (Color online) Determination of the drift between two successive mean flow patterns 
in order to quantify the slow variation, (a): Magnitude of the difference between two successive 
mean flow patterns taken at zq = H/8 with r = 50Tf. (b-d): Vfe+i,fe(2:o) versus averaging interval 
k taken at three different zo which are indicated in the legend and the same for all three panels. 
Data are for run 2 in table(b) r = lOOT/, (c) r = SOT/, (d) r = SOT/. 


that is by two orders of magnitude smaller. A time average taken over t has to be long 
enough such that the turbulent fluctuations in the velocity field are suppressed {t ^ Tf). 
However, if the averaging procedure proceeds over a very long time interval then these 
patterns will be washed out for all Rayleigh numbers discussed here. We can decompose 
the velocity and temperature fields into a time-averaged field and remaining turbulent 
fluctuations as 


Ui{xj,t) = {ui{Xj))tu'i{xj,t) and T{xj,t) = {T{xj))t-\-T'{xj,t) , (3.1) 

where {■)t denotes a time average. In Fig. we analyze the slow drift of the large-scale 
flow pattern. The total integration time interval, T = Mr, is divided into M equidistant 
subintervals Ik with k = 0, M— 1. These averaging intervals are taken from kr to (fc-l-l)r 
with T Tf. Panel (a) of the figure shows the magnitude of the difference between two 
successive mean flow patterns taken dX zq = H/8. We see that pointwise differences get 
as high as O.bUj in this example. Panels (b)-(d) of Fig. [^show a measure for the drift 
of the mean flow patterns which is defined as 




Mzo))- Mzo))A.tei, 


(3.2) 


The notation {■)A,t stands for a plane-time average. Run 2 is advanced for 400T/ and 
this time interval is split into fractions of r = 20,50 and lOOTf, respectively. While for 
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Figure 3. Streamline plots of the velocity field in the Rayleigh-Benard convection cell at 
Ra = 500 000 and Pr = 10. A view from the top onto a quarter of the cell is displayed. Left: 
instantaneous streamline snapshot. Right: time-averaged streamline plot obtained for averaging 
over T — 50Tf. Data are from run 5 in table 


T = lOOT/ the drift velocities in all three planes are of the same size, the data for the 
two smaller r imply that the drift in the center plane is in parts slightly slower. For 
averaging times r ;< 20Ty we reach the range of typical turnover times of a Lagrangian 
tracer within a large-scale circulation roll (Emran & Schumacher (20101). Therefore, we 
do not consider smaller time intervals r. In all three plots, we detect nearly the same 
magnitude of the drift velocity. This allows us to derive a time scale of the processes 
which is H/Vk+i^k ^ lO^T/, a large time scale which is not accessible in this study. 
This time scale is comparable to that of a slow spanwise drift of streaky structures in 
plane Poiseuille flow which has been reported very recently by Kreilos et al. (2014). This 
estimate is also consistent with the one for a time scale of horizontal motion, T/j, that 
should vary as = T^Tj. Furthermore, we observe that the drift velocities for planes at 
z = St, H/2 and H — St are of same order of magnitude. This suggests that the mean 
flow roll pattern drifts slowly as a whole. 

Figure repeats the analysis at Pr = 10 and Ra = 500 000. Now, the streamlines of 
the snapshot appear much less disordered than for Pr = 0.7. Consequently, the difference 
to the mean flow pattern is much smaller. One reason could be that the thermal diffusion 
is less compared to the momentum diffusion when the Prandtl number grows for a fixed 
Rayleigh number. This results in thermal plumes which have thinner stems and disperse 
less rapidly with respect to time. Thus the stirring of the fluid by plumes is less efficient. 
The result is in line with the decrease of the Reynolds number for growing Prandtl 
number as shown in Tab. Our finding is also supported by Silano et al. ( 2010| who 
have observed decreasing peak velocities for increasing Pr. 

In Fig.|^ we summarize the results of the Reynolds-decomposed velocity field (see the 
decomposition in Eq. (3.11). In detail, we define 


^rms — 




(3.3) 


We include further runs at the same resolution which are not listed in Table but 
in the caption. The smallest Rayleigh number was Ra = 2 000 for Pr = 0.7 which is 
slightly larger than the linear instability threshold, Rac = 1 708. When expressed as a 
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Figure 4. (Color online) Root mean square values of the total velocity, the time-aver¬ 
aged velocity a nd t he remaining turbulent fluctuations as a function of Rwleigh number 
Ra (see Eqns. (3.3|). Additional data points beside those listed in Tab. [u are given at 
Ra = 2 000,3 000,4 000,10 000,20 000,50 000 and 100 000 for the series at Rr = 0.7 and 
Ra = 50 000 for Pr = 10. 


distance to the linear instability threshold this gives e = {Ra — Rac)/Rac = 0.17. In 
this case, velocity fluctuations are practically absent, the flow pattern is almost steady 
consisting of several subdomains with stripe textures. With increasing Rayleigh number 
fluctuations of all three parts of the velocity field (see Eqns. (3.3)) grow up to Ra « 5 000 
which corresponds to e = {Ra — Rac)/Rac = 1.93. At about this Rayleigh number, Urms 
reaches a local maximum and starts to decrease with increasing Rayleigh number. At 
Ra « 10 000, the turbulent fluctuations Vrms exceed Urms- At Ra ^ 100 000, Urms and 
Vrms reach a local maximum and level off. For this Rayleigh number, the flow is already 
turbulent, the fluctuations Vrms are by a factor of two larger than Urms- We also show 
three data sets for the case of Pr = 10. The magnitudes of all three parts are significantly 
reduced which confirms our observation from Fig. Up to the accessible Ra = 500 000 
all three terms continue to grow suggesting that the maxima are shifted to higher Ra. 


3.2. Estimate of turbulent viscosity and diffusivity 

The next step is to estimate the turbulent viscosities and diffusivities in the bulk of the 
cell and to evaluate the resulting turbulent Rayleigh and Prandtl numbers. We start with 
the Boussinesq ansatz for the closure which connects turbulent fluxes (or stresses) with 
the mean gradients (see e.g. Wilcox (2006); Shams et al. (2014|)) and states that 


dui 


u'u! = 

* dx, 


u'T' = -kI^ 


, dT 


dx-i 


(3.4) 


where bars denote an appropriate space-time average. Our following estimate will aim 
at obtaining numbers u* and k* rather than exploring the full tensorial structure of the 
turbulent viscosities and diffusivities. This would go beyond the scope of this work. We 
will restrict the analysis to the dominant contributions only. 

In case of the turbulent diffusivity, we focus to the vertical transport of heat from 
the hot bottom plate to the cold top plate. The comparison of the three convective 
fluxes shows that the magnitude of the mean vertical flux is the largest. The turbulent 
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Figure 5. (Color online) Vertical profiles of the plane and time averaged correlations and 
derivatives which are required to determine the turbulent viscosity and diffusivity. Data displayed 
in the figure are obtained for Ra = 500 000 and Pr = 0.7. 


diffusivity, k*, can be obtained by the following bulk average 

/uv/ d{T{z))A,t , 

/ KuJ^{z))A,tdz = -n^l -^- dz. (3.5) 

Jst Jst ^z 

Time-plane averages are denoted by {■)A,t- Figure [^displays the resulting profiles which 
enter the determination of k* via (3.51. The double-headed arrow indicates the bulk 
region in the figure. 

In case of the momentum transport the determination is less straightforward. We 
can expect that the horizontal turbulent mixing is also important. First, we proceed 
however similar to the temperature field. We take the magnitude of horizontal velocity 
uj_(r, ((), z, t) = Ur{r, (p, z,t)er + p, z,t)e^. This field is decomposed again into a 
temporal mean and remaining fluctuations. The turbulent viscosity, u* is determined in 
a similar way as the turbulent diffusvity 

d{u^{z))A,t , 

/ {u^Uj_(z))A,tdz =-ly^ --- dz. (3.6) 

J S'T 'y C/Z 


Here u'^ and u±_ denote magnitudes. The resulting profiles that enter ( |3.5[ ) and (3.6) are 
displayed in the right column of Fig. In Tab. [^we summarize the resulting turbulent 
Rayleigh and Prandtl numbers which result from this closure procedure. The turbulent 
Rayleigh numbers, i?a*, are reduced for all three cases. Ra* gets consistently smaller 
with decreasing Prandtl number Pr since the amplitude of the turbulent fluctuations 
increases. We obtain Pr^ < Pr for all three cases. Their magnitudes vary between 0.2 
and 0.4. The turbulent Prandtl number Pr* increases slightly with increasing Pr. 

In case of runs 2 and 5, we then conducted a DNS with the same molecular viscosity 
and diffusivity as i^* and k*, respectively. The resulting streamline pattern for run 2 is 
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Run 

Ra 

Pr 

Rat 

Prt 

^0 


b*o 

bl 

2 

500000 

0.7 

4 500 

0.21 

15 

8 

18 

7 

5 

500000 

10 

39 000 

0.38 

20 

1 

18 

1 


Table 2. Betti numbers bo and bi for the original simulations at Ra and Pr as well as 6 q and 
6* for the corresponding runs at i?a* and Pr*. The number of the runs corresponds with Tab. 
Temperature patterns at mid plane have been analyzed. 



Figure 6. (Color online) Streamline plots of the velocity field in a Rayleigh-Benard convection 
cell as a view from the top. Left: plot averaged for 200 T/. Right: magnification of the same 
data. The DNS was conducted for Ra — 4 500, Pr = 0.2 which corresponds with the turbulent 
viscosities and diffusivities that correspond to the time-averaged data in Fig. [i](&h). 


displayed in Fig. A time average over 200 Tf was applied at Pr = 0.2 and Ra = 4 500. 
The flow structure has to be compared now with the time-averaged one from Fig. 
(g,h) and indeed a reasonable visual agreement of both large-scale patterns is found. 
We determined the Betti numbers {6o,5i} from two-dimensional horizontal cuts of the 
mean temperature at z = 1/2 in both cases (Kurtuldu et al. (20111). Betti numbers 
are d positive integers to characterize a d-dimensional set topologically. In detail, bo is 
the number of connected filaments which is obtained by digitizing a grayscale picture 
at a threshold, bi counts the number of enclosed holes in the pattern. We choose the 
temperature field in the mid plane. A threshold temperature T = 0.5 results in Betti 
number pairs which are listed in Tab.j^for runs 2 and 5 as well as their corresponding runs 
at i?a* and Pr*. Additionally, we estimated the average width of the rolls by counting 
the mean number of rolls that fit into the cell along different orientations. The values 
vary always around a width of 2H, but are not exactly equal. 

We also determined the turbulent viscosities from horizontal turbulent diffusion pro¬ 
cesses. It turns out that a simple adaption of the averaging procedure of Eqns. (3.5) and 
(3.6) to a radial dependence is not successful. The roll patterns cause radially oscillating 
profiles which result in strong cancellations for the averaged turbulent stresses and mean 
strain rates. If we omit the radial averaging and analyze the local Boussinesq relation 
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{u'^u'j{z))^^H- 2 ST,t = -v*{'r)d{uj{r)) ^^h-25t, tl dr for j = r, (j),z, we get indeed turbulent 
Prandtl numbers which are locally closer to one, but vary significantly with r. 


The resulting i?a* and Pr* are such that the DNS yield time-dependent patterns, 
in particular for run 5 with Pa* = 39 000 and Pr* = 0.4. We therefore repeated this 
“renormalization procedure” in the weakly nonlinear regime and obtain Pa** = 3000 
and Pr** = 0.13 for run 2 and Pa** = 1800 and Pr** = 0.17 for run 5, respectively. 
Both runs end thus in the convection regime close to the onset. 


To summarize this section, all routes of analysis will in general not lead to turbulent 
Prandtl numbers Pr* « 1. A turbulent Prandtl number smaller than unity can be in¬ 
terpreted as follows: plume filaments of the temperature are coarser and diffuse faster 
than vortex filaments next to them. This circumstance could be connected to the fact 
that the width of rising and falling plumes is of the size of the thermal boundary layer 
thickness which is rather large for our Pa. In contrast, vorticity is frequently generated 
on finer scales. Vortex filaments are for example generated by locally reversed flows next 
to rising plumes, a consequence of incompressibility. An increase of the Rayleigh number 
to very large values could then increase the turbulent Prandtl number to one since the 
typical flow structures are getting finer and the boundary layers themselves are expected 
to become eventually fully turbulent. 


At this point, it should also be mentioned that the particular magnitudes of turbu¬ 
lent Prandtl numbers, Pr* are still an open problem. For example, Spiegel (19711, Kays 


(1994), or Grotzbach (2011) discuss the dependence of Pr* on the distance from walls or 
on the original Pr. In case of homogeneous isotropic turbulence, Nakano et al. (19791 de¬ 
rived a value of Pr^ = 0.4 from a spectral formulation based on the classical Kolmogorov 
turbulence theory. 


3.3. Robustness of large-scale mean flow patterns to additional side wall forcing 


The sensitivity of SDC patterns to side wall effects and suppressed mean flows has been 
discussed in Bodenschatz et al. (1991) and Chiam et a/. | ( [2003 ), respectively. This moti¬ 
vates us here, also in view to spontaneous symmetry breaking, to study their robustness 
with respect to an addition of a volume forcing to (13. The forcing is set up such that 
it sustains a steady Lamb-Oseen-type vortex Ui{r, z) with i = {r, z} very close to the 
side walls of the cell at (rp = (P — l)/2, zq = 1/2). This vortex generates an azimuthally 
symmetric mean flow at the side wall. The circulation, n, and the radius of the vortex 
core, tl, are chosen such that no-slip boundary conditions can still be satisfied by setting 
this flow to zero below a certain threshold. This clearly prohibits a stronger variation of 
the amplitudes and thus of the strength of the additional forcing. 


Incompressibility of the full velocity field is sustained via the solution of the Poisson 
problem for the pressure in each time step. Figure shows the results for the mean flow 
pattern. In case of Pr = 0.7, the toroidal roll is clearly visible right at the side wall. 
A second roll next to the side walls can be established by the additional forcing term. 
Towards the center of the convection cell the mean pattern remains however unchanged 
as can be seen by a comparison with panels (g,h) of Fig. The turbulent fluctuations 
are large enough to re-establish the mean flow pattern. This is different for Pr = 10. In 
comparison to Fig. the pattern has changed significantly. The toroidal roll pattern of 
the time averaged velocity is continued almost to the center of the cell. The reason for 
the stronger impact of the additional side wall forcing lies in the significantly lower level 
of turbulent fluctuations which we documented in Fig. 
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Figure 7. Streamline plots of the velocity field in a Rayleigh-Benard con vect ion cell as a view 
from the top for Ra = 500 000 with the additional forcing fi (see also Eq. (2.11). Left: Pr = 0.7. 
Right: Pr = 10. Both data sets have been averaged over 150Tf. We took il — 1 for the (non-di¬ 
mensional) circulation and tl ~ 0.01 for the (non-dimensional) radius of the vortex core. 


4. Summary 

We presented three-dimensional DNS of thermal convection in the soft turbulence 
regime to study time-averaged velocity field patterns and their dependence on the Prandtl 
number in very large aspect ratio convection cells. Our DNS demonstrate clearly that 
the SDC patterns, which are known from the weakly nonlinear regime, continue to exist 
in the turbulent regime. They remain thus dynamically relevant and do not simply dis¬ 
appear when convection turns into the turbulent regime. The patterns are revealed when 
the turbulent fluctuations are removed by time averaging over intervals of the order of 
lO^Tf, which is significantly smaller than the time scale over which the mean velocity 
and temperature patterns evolve. Our simulations allow us to calculate the turbulent 
viscosities and diffusivities as well as related turbulent Rayleigh and Prandtl numbers, 
and Pr*. Their values fall indeed back into the range of the original SDC regime. 
The turbulent Prandtl numbers Pr* vary between 0.2 and 0.4 and increase with increas¬ 
ing Pr. Our studies showed also that the mean patterns are robust to finite-amplitude 
perturbations once the turbulent fluctuations in the flow are sufficiently large, i.e., once 
Pr at a given Ra is sufficiently small. We demonstrated this by a side wall forcing that 
sustained an azimuthally symmetric vortex. 

Three future implications follow to our view: (i) it has to be investigated systemati¬ 
cally if the mean flow patterns which are similar to SDC persist to even higher Rayleigh 
numbers or if the mean flow structure is changed. This would require numerical studies 
at high Rayleigh numbers and large aspect ratios, (ii) a more detailed analysis of the 
turbulent viscosities and diffusivities for larger Rayleigh numbers will provide useful in¬ 
put for technological and astrophysical applications in which the small-scale convective 
turbulence has to be modeled. This would however imply to explore systematically the 
tensorial nature of the turbulent viscosity which we did not analyze in the present work, 
(iii) our results could also provide useful input to reduce the degrees of freedom system¬ 
atically and to derive some effective equations for the large-scale patterns, as done in 
other systems (Malecha et al. (2014[)). 
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